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Abstract 

We demonstrate robust and reliable signatures for the transition from quantum to classical be- 
havior in the position probability distribution of a damped double-well system using the Qunatum 
State Diffusion approach to open quantum systems. We argue that these signatures are within 
experimental reach, for example in a doubly-clamped nanomechanical beam. 
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I. INTRODUCTION 



The difference between quantum and classical behavior is of fundamental interest with 
many experimental consequences in contexts such as superconducting quantum interference 
devices (SQUIDs), cold atom systems, and significantly in nanoelectromechanical systems 
(NEMS). NEMS are small, oscillate at high frequencies and are typically maintained at low 
temperatures. Most of their parameters are adjustable, making them an ideal setting to 
explore the transition from quantum to classical behavior. Recent experimental results 11, 
12J, indicate that such systems can now be prepared and observed in their ground state, 
thus enabling the study of genuinely non-classical behavior in the dynamics. The exciting 
prospect of being able to observe the quantum-classical transition in detail has contributed 
to interest in quantum mechanical nanoelectromechanical systems (occasionally refered to as 
QEMS in the literature), whose study is an active and rapidly expanding field of research^. 

A seminal theoretical study of the quantum-classical transition in NEMS was carried 
out by Katz et al|4| using a nonlinear resonator, which is a paradigmatic model for several 
common types of NEMS that yields quantum behavior distinct from the classical. They 
considered an isolated resonator as well as one modeled as an open system (coupled to the 
environment) in search of experimental signatures of such quantum behavior. Specifically, 
they examined the dynamics of the classical phase space distribution in such a system com- 
pared to that of the quantum-mechanical Wigner function, and looked at a zero temperature 
and a finite temperature environment. They were able to find differences in experimentally 
accessible signatures: In particular, the quantum mechanical version of the resonator has 
non-zero probability of being found in a position where the classical resonator has zero prob- 
ability of being found. However, these differences are relatively small and blur away quickly 
when experimental noise or finite temperature effects through thermal noise are considered. 

In this paper, by considering a slightly different but equally accessible nonlinear oscillator, 
we are able to obtain more pronounced experimental signatures of the quantum-classical 
transition; these are arguably more visible at finite temperatures and in the presence of 
experimental noise. Further, we show that rather than comparing just quantum and classical 
behavior, there is insight to be gained from studying the continuous transition where the 
classical behavior emerges as the limiting case of the quantum dynamics. That is, although 
systems are fundamentally described by quantum mechanics, their behavior changes as we 
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increase the size of the system (affecting the characteristic action scaled by %), the effect 
of the environment (decoherence) or other parameters of the system. For our system, we 
indeed recover the classical behavior as the limiting case. The details of the emergence of 
the classical limit are illuminating, and in particular, the signatures do not always change 
monotonically, thus establishing that comparing a single quantum solution with a classical 
solution leads to an incomplete picture. 

Katz et al considered the sytem with the Hamiltonian 

H S vs = ~P 2 H — x 2 H — x 4 — xFcosut. (1) 
sys 2 2 4 

This quadratic nonlinear oscillator is within experimental reach in the near-quantum regime, 
and is an excellent system to study. However, even more interesting effects such as quantum 
tunneling obtain if we change the sign so that the x 2 term is negative^. This yields the 
so-called double- well Duffing oscillator, which has the Hamiltonian H sys = ^p 2 — |x 2 + 
— xF cos cut, and allows for clearer signatures of the classical to quantum transition as 
we demonstrate below. 

This paper is organized as follows: We start by motivating the modeling of a specific 
NEMS as a doubly-clamped beam to show the connection between experimental paramters 
and model parameters. We then sketch the quantum state diffusion (QSD) approach, where 
Lindblad operators act within stochastic Schrodinger equations to incorporate the effect of 
the environment, used to model the behavior of such an oscillator understood as an open 
quantum system. These two allow us to establish that changing the sign of the quadratic 
term is within current experimental capabilities and that such experiments are in or close 
to the realm where the quantum to classical transition can be explored. We then present 
results and conclude with our analysis. 



II. NEMS MODELED AS A DOUBLY-CLAMPED BEAM 



Many nanoelectromechanical devices can be modeled as doubly-clamped beams driven 
near resonance and continuum mechanics continues to serve as an adequate model even at 
the submicron scale 0-3. Usually magnetomotive [9] or optical [6] actuation is used to study 
the resonance behaviors of NEMS such as dynamically induced bistability, hysteresis and 



effects of parameter attenuation on energy dissipation. The resonant behavior of the doubly- 
clamped NEMS structure has been experimentally identified with the fundamental bending 
mode [h]] of the problem. We sketch here the theory allowing us to relate the resonance 
frequency to characteristic parameters of NEMS and to determine the scale of the system 
with respect to H, establishing the regime of the transition to classical behavior in these 
systems. 

Starting with the assumption of an ideal beam of rectangular cross-sectional area, we can 
obtain the potential energy V of a do ubly -clamped beam under tension along the longitudinal 
direction. This consists of two parts 11] : elastic and bending potential energy, 

Veiastic = T([° v 7 ! + (V') 2 dx - l ) (2) 
Jo 

and 

EI r l ° a 

Vbending = — / (y'fdx. (3) 
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T is the external force on the beam in the form of the tension in the longitudinal direction 
along the beam, E is the elastic modulus, / is moment of inertia of the cross-section and Iq 
is separation of two ends that clamp the beam. Note that Iq need not be the length of the 
beam when it is neither stretched nor compressed because a static doubly-clamped beam 
may already be buckled due to inherent tension along the beam. V e i as uc is due to tension 
built up in the beam, which causes the elastic beam to deform and therefore gain energy. 



This is usually referred to as strain energy 



121 ] while Vb en ding is purely due to bending. The 



quantity EI, termed "flexural rigidity", denotes the force necessary to bend a beam by a 
unit of curvature. 



The tension T along the beam has two parts 13|] : firstly, inherent tension To (either from 
compression at both ends of the beam or manufacturing process) and additional tension AT 
(due to transverse motion that stretches the beam), whence T = To + AT. When the beam 
is stretched, its length changes and therefore the additional tension AT is given by fractional 
change in length multiplied by elastic modulus and cross section area A|l3j:AT = l -^-EA, 
where I refers to the actual length of the buckled beam. We can rewrite this length as 

l = l + Al= f° dx^Jl + {y'Y « f° dx(l + \{y'f) = Iq + \ / ' \y'fdx. (4) 
Jo Jo * * Jo 



This yields that the change in tension AT = ^ J^°(y'fdx whence the total tension 

f A r l ° 

T = T + — ^ (y'fdx. (5) 

Substituting these expressions for Al and T in the expression for elastic potential energy 
above, we get 



T f l ° FA f l ° 

V elastlc = -f I (y'fdx + =£[] (y'fdx] 2 . (C)) 



Thus, the total potential energy is 

v = v elastic + v bending 

= f ti\y>fdx + ^(y'fdxf + f tfiy'Tdx. 

At low temperatures, as is often true for NEMS experiments, it is sufficient to consider only 
the fundamental mode. Mode amplitude solutions for this fundamental mode y(x) are of 



the form 14( y(x) = \Y\l — cos(27Tp)], where Y is the transverse displacement of the center 
point of the beam. When we make this substitution, the fundamental mode amplitude Y - 
that is, the location of the vibrating central point - is effectively the position variable of a 
single particle evolving in a potential V, which as a function of Y is given by 

v ( y ) - -i^ Y +^r y (8 ) 



i EAtj 4 
16/3 • 



Here /, the moment of inertia of rectangular cross section is ^ 15], where a is the thickness 



of the beam and b is its width. I varies with the geometry of the beam for specific experi- 
mental setups. Further, depending on whether inherent tension T compresses or stretches 



the beam, the sign of T can be negative or positive. Compressive inherent tension a 



ong 

3- 



the beam can arise from manufacturing processes, and can vary from 200 MPa to 1 GPa 
In general, the inherent tension is not known a priori in an experimental setup. Instead, by 
considering the continuum mechanics model and fitting the experimentally measurable res- 



onance frequency, the inherent tension can be calculated [lOJ]. For sufficiently large inherent 
compressive tension, T < — 47r !j E J , we obtain a double- well potential and the beam is in a 



state of "Euler instability" |8|, [16|], meaning that it is buckled. 
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Thus a single-well potential or a double-well potential can be experimentally achieved in 
doubly-clamped NEMS structures and these systems can then be studied for their dyamical 
behavior. The dynamics of a single- well NEMS actuated near the fundamental bending mode 



frequency, has been extensively studied 



yet been fully explored. The results of Ref. 



18 1 but that of the double- well case has not been 



10[ suggest that it within experimental reach 



and the above clarifies that in general tuning the characteristic parameters should enable 
switching from a single-well to a double-well potential. 



III. OPEN QUANTUM SYSTEMS 

We now turn to the quantum version of the model for the problem and connect char- 
acteristic parameters with the degree of 'quantumness' of the system. An open quantum 
system with weak system-environment coupling and a Markovian environment is modeled 
by a master equation 



m =— [h, p] - \ + p^h 
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(9) 



where Lj is the Lindblad operator representing the system-environment interaction and p 
is the reduced density operator. The quantum state diffusion (QSD) approach corresponds 
to solving quantum trajectories which are unravellings of the master equation. This al- 
lows for numerical efficiencies, and physical insights not available via the master equation. 
Specifically, we can use the QSD numerical library 
Schrodinger equation 



19] to solve a stochastic version of the 



d\if)) =^H\i))dt 



+ E((4V J -^4%-^(4 t )(4)) 



\ip)dt ^ 



The solution \ip(t)) to Eq. ffTUl) is a "quantum trajectory", and we obtain p = jj 
as the mean over an ensemble of M normalized pure state projection operators 19J . This 
p satisfies the Lindblad master equation Eq. fl9]). That is, starting with an ensemble of 
identical pure states, we obtain the time evolution of multiple trajectories that evolve into 
an ensemble of different pure states due to interaction with the environment; this is the 
density matrix. 



For t 



as 



20 



re open quantum double-well Duffing oscillator H and Lj in Eq. f jTOj) are chosen 
22 1 



H = H D + H R + H e 



muj. 



i 

H R = ^(xp + px) 



— V + 

7, 



^4 



x — 



°x 2 



- — gmlulx cos (cut) 



mu 



-x 
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rP 



(11) 
(12) 

(13) 

(14) 

(15) 

(16) 



where 7 is the strength of coupling between the oscillator and thermal bath environment. We 



now define x and p in unitless forms Q = \J muo / hx and P = ^Jl/moJohp respectively 2ll.l22] 
yielding the dimensionless set of equations 



H' 



D 



H 'ex 

L'i = 
L' = 



-I^QP + PQ 

■- — j;Q cos(Oi) 



/T(l + n) [Q + iP 

/ rn(Q-iP 



(17) 

(18) 
(19) 

(20) 
(21) 



where the time in measured in units of uj^ 1 so that T = •y/ujQ,Q = u/loq. In Eqs. (|15[ 
[T6l20|2ip . n = (e n ^/ k BT — \\ 1 is the Bose- Einstein distribution of thermal photons rep- 
resenting the effect of the environment at finite T and is evaluated at the natural fre- 



quency of the oscillator 



23] . This dimensionless formulation isolates the scaling factor 



[3 = y ^ as the single length scale of the problem, where S is the classical action, such 
that (3 2 = m ^ WQ ■ For the classical equation of motion for a damped Duffing oscillator, which 
* s It? + ^lt ft 1 ^ ~ x = § cos (^); solutions to the classical equation of motion are un- 
changed with respect to different values of (3 except for the length scale of the phase space. 
However, changing (3 changes the quantum solution considerably. 

Comparing the quartic potential Eq. (jSJ) for the motion of the fundamental mode of the 
nanomechanical resonator with the potential term in the Duffing Hamiltonian Eq. (1121) . it 
is straightforward to see that for the quantum nanoelectromechanical resonator 



2 _ fi hlo A 

= s = s^y 2£p/3 | ( _ 1)(1 + _zj|_ )]S 

where, as noted in section II, we restrict our consideration to resonators of the double well 
shape with T < — ^■F 1 ■ The scaling parameter (3 2 determines the way in which the action 
of the nanomechanical system scales relative to h. 

As such, studying the system as (3 is changed allows us to explore the quantum to classical 



transition. Previous theoretical work 
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^22] for the Duffing problem focused on the behavior 
of the quantum trajectories has confirmed that — > indeed recovers the classical limiting 
behavior, including the presence of 'strange attractors' of chaos in the phase space, while 
(3 = 1 shows behavior distinctive of the quantum regime. There is interesting physics in 
the intermediate regime (20], and in particular the details of the transition as f3 changes 
are informative (as below). As shown in Table (jl]), doubly-clamped nanoelectromechanical 
beams can be , sinele-walled carbon nanotubes 



. The analysis 
2511 . which has 



(SWNT) or multi- walled (MWNT) nanotubes |6| and even metals like gold|24f 
above for the Duffing oscillator also works for a doubly-clamped Pt nanowire 
a high fundamental frequency (greater than 1 GHz). The characteristic parameters given in 
Table (jl]) for example show that the current experimental setups can be tuned within the 
range of (3 considered in this paper. Specifically, for the SWNT experiment of Witkamp et 



allOj], (see Table HJ, they used a device of radius r = l.Qnm and length Iq = 1.15/im. This 



4 

yields, using their other device parameters and an I = ^ that 



(3 2 = 5.34X1(T 6 X(-^- - l)- 3/2 (23) 



where A = -^f- is greater than Air 2 to yield the double-well shape. If A is in the range 
(39.5 — 60), we get /3 in the range (0.65 — 0.04), which as we see below, is precisely the range 
needed to map the transition. Witkamp et al report A = 26, so this is clearly well within 
reach of current experiments. 



TABLE I: Four common types of doubly-clamped structure of NEMS that have been 
realized experimentally. The characteristic parameters can be used to find the fundamental 

mode frequency. 



Type 


Length, / ( nm ) 


Diameter (nm) 


Density p (kg/m 3 ) 


Elastic ModulmE(TPa) 


Si[5j 


50-10 4 


5x10 to 10x20 


2330 


0.137 


SWNT[12],[g] 


50-10 4 


1 


1930 


1.25 


MWXT[6J 


50-10 4 


20 


1930 


1.25 


Pt nanowire[25l 


1300 


43 


168 


21090 



There are other parameters of interest related to the system-environment interactions. 
In particular, the experimentally measurable quality factor Q quantifies energy lost in com- 
parison to total energy of the resonator during one complete driving period which can be 
related to the damping T. A comprehensive account for the mechanisms of quality factor, 
Q, in NEMS has not been established 24|. but for the case of doubly-clamped beam, Q 
can arise from clamping, thermal elastic damping, as well as defects of crystalline struc- 
ture of the material [2| . Usually, Q can be measured experimentally using magnetomotive 
techniques j^]. For low temperatures (20 K and below) Q can be as large as of the order of 
10 3 [l3]. Lowering the temperature further to the milikelvin range further depresses energy 



dissipation and recent experiments 24j, |26J| have achieved quality factors on the order of 10 5 . 



Although these experiments have pushed towards high Q and correspondingly low T, it is 
clearly possible to make systems more dissipative, thus allowing the exploration of a large 
range in T. 



IV. RESULTS AND DISCUSSIONS 



The solutions to Eq. ( TlOj) yield "quantum trajectories" which have been previously 
discussed in detail in the context of the recovery of classical behavior from quantum 



mechanics 



20H22j . What we focus our attention on below is ensemble-and-time-averaged 



asymptotic behavior of the system, defined through the probability distributions 

M N 

=EE^* (t)l ^ )> " (24) 
i=l k=l 

with ^transient + 47r > t > transient- in doing this we are averaging over multiple M tra- 
jectories generated in QSD simulations to yield the behavior of density matrices accessible 
the laboratory. We also eliminate phase-dependent idiosyncracies of these distributions and 
correlate better with experimental capabilities regarding the measurement of the dynamics 
of the quadratures by time-averaging trajectories \ipl(t)) with N samples of each trajectory 
over two driving periods, taken after a suitably long transient time transient- 

The time averaging is not necessary but is useful because (a) in exploring the quantum to 
classical transition, the properties that signal change most clearly as parameters are varied 
are asymptotic, global characteristics of the system rather than unrepresentative "snap 
shots" at specific times. Also (b), the time-averaging eliminates experimental difficulties in 
determining the phase of the driving relative to the observation, as well as allowing one to 
quickly build statistics. 

We have explored P aV g{x) at both zero temperature (that is, n — 0) and finite temperature 
for several /3; we report here /3 = 0.01,0.3 and 1.0. For each of these values, as we change 
the coupling strength T, we also looked for signatures of the quantum to classical transition 
through the changes in P aV g(x). Fig. (JTJ shows P a vg(x) at zero temperature. With high 
dissipation (r = 0.3) we see what we term nonmonotonicity in the behavior of the probability 
distribution as a function of j3. Specifically, the probability distribution changes as follows: 
(a) At /3 = 0.01 it is an asymmetric object centered in a single well, followed by (b) at 
(3 = 0.3 an essentially symmetric double-peak structure across both wells and then (c) at 
P — 0.1 a symmetric single peak straddling the central potential barrier. 

We understand these results as follows: For (3 = 0.01, the most classical case, the dissipa- 
tion results in the particle being unable to overcome the central potential barrier even with 
driving and therefore being confined in a single well. The probability peaks at both ends of 
the probability distribution correspond to the classical turning points of the oscillator orbit 
in phase space. That is, the dynamics in this situation are identical to those that would 
be generated by a classical trajectory, and unsurprisingly this means that the probability 
distribution is also entirely explained by classical considerations. 
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As we increase the quantum-ness of the system (decrease the scale of the system) to 
(3 = 0.3, the particle acquires a nonzero probability of being in both wells. This is classically 
impossible given that it is a zero-temperature problem with high damping. The inter-well 
transitions could indeed occur because, as the system gets more quantum-mechanical with 
(3 increasing, the scale of the Lindblad term effectively also increases. It is therefore possible 
that the particle continues to act essentially classically but the increased effect of the noise 
leads to a classical noise-activated hopping over the barrier. However, our studies show 
that the particle is behaving entirely quantum-mechanically during the interwell transition, 
whence this transition is a quantum tunneling effect. Specifically, we have seen that the 
Wigner function for the particle shows classically forbidden negative regions throughout its 
dynamics. We have also seen that the energy expectation value remains negative while the 
expectation value of the position transits past the origin 20( . Finally, we have found instances 



in time when a wave-packet has significant probability on either side of the barrier, which is 
incomensurate with it being a localized wavepacket behaving like a classical particle hopping 
over the barrier. Incidentally, the last two demonstrate one advantage of using the QSD 
wavefunction-based approach, since such information is unavailable in the master equation 
approach. The symmetric double-peak structure of P aV g(%) at this value of (3 indicates 
that the ensemble average over the wave packet behavior leads to an asymptotically equal 
probability in either well. 

For (3 = 1, the most quantum case, the probability is peaked at the origin. There are 
multiple ways of understanding this. The first is that the zero point energy is such that the 
lowest states of the potential are above the barrier. As such the probability is peaked in the 
middle of the quartic larger well, ignoring the quadratic bump at the bottom. Another way 
of understanding this situation is to realize that as (3 is increased above 0.3, the scale of the 
system changes such that the two peaks of the symmetric distribution from (3 = 0.3 come 
closer together until they merge. 

These three figures show not only that there is a clearer difference between the quantum 
and the classical probability behavior in the double-well Duffing oscillator compared to the 
single-well version, there are interesting things to be learned from studying the intermediate 
regime. 

As as alternative case, we also present results from slightly lower dissipation, with T = 
0.125. In this case, even at the classical limit, the dissipation is not strong enough to confine 

11 



the particle in one wel 
across both wells 



. Specifically, the particle trajectory is that of a strange attractor 
20M22| which yields the unusual shape for the probability distribution. The 
difference in the classical behavior does not persist, and the distributions at higher for 
this T value follow the same pattern as those for V = 0.3. As such, comparing the two sets 
of figures, the results at V — 0.3 show a more dramatic transition from the classical to the 
quantum regime. 



i,r, n=o 



r=0.125 



r=o.3 



£=0.01 




V 



8=0.3 





(8=1 





FIG. 1: P avg (x) plots for various parameters. The degree of quantumness increases from 
left to right as 0.01, 0.3, 1.0. The values for T = 0.125, 0.3. 



At even lower values of T, the vanishing dissipation means that the wavefunctions do 
not localize in a basis as happens for greater dissipation 19( . This leads to the standard 
semiclassical convergence problem for the calculations of P avg (x). However we are able to 
compute P a vg(x) at finite temperature for each set of (3 and T and these are presented 
in Fig. ((2]). We see quickly that the results are similar to those at zero temperature in 
Fig. fll]) though we can see that temperatures corresponding to ~ 4 thermal photons begins 
to eliminate the structure of the distribution(s) at (3 = 0.3; this is roughly where quantum 
effects are equal to thermal effects. As we can see, for large dissipation (r = 0.3) and the 
most classical case, the particle is still confined in one well; that is, thermal activation is not 
high enough to induce switching between wells. The probability distributions are of course 
broadened due to thermal effects, but the signature of the transition as /3 is increased are 
robust. The last row of Fig. fl2]) reassuringly shows the existence of a smooth transition 
away from zero temperature behavior. 

In conclusion, we have shown that in a open double-well quantum oscillator clear sig- 
natures of quantum behavior and in particular the transition away from classical behavior 
can be found in the changing shapes of an asymptotically obtained probability distribution 
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FIG. 2: P avg {x) plots for various parameters at finite temperature. For the first three rows, 
the number of thermal photons n = 4 and /3 increases from left to right as 0.01, 0.3, 1.0. T 
increases down the grid as 0.03, 0.125, 0.3. The last row has T = 0.3 and n = 0.5 and shows 
a smooth transition from the zero temperature case in Fig. (JTJ) . 

P aV g{x)- At the right parameters, this transition from classical through quantum mechanical 
behavior can be nonmonotonic and hence should be clearly visible. Further, these transi- 
tions are robust at finite temperatures. The experimentally tunable parameters in NEMS 
allow for flexibility in exploring the parameter landscape of the transition and are excel- 
lent candidates to study this transition. We believe that such experiments are viable in 
doubly-clamped NEMS, (see Table]!]) 
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